arXiv:1501.01069vl [astro-ph.HE] 6 Jan 2015 


Astronomy & Astrophysics manuscript no. paper' cheers' v2 © ESO 2015 

January 7, 2015 


Chemical Enrichment RGS cluster sample (CHEERS): 
Constraints on turbulence 

Ciro Pinto 1 , Jeremy S. Sanders 2 , Norbert Werner 3,4 , Jelle de Plaa 5 , Andrew C. Fabian 1 , Yu-Ying Zhang 6 , 
Jelle S. Kaastra 5 , Alexis Finoguenov 7 , and Jussi Ahoranta 7 

1 Institute of Astronomy, Madingley Road, CB3 OHA Cambridge, United Kingdom, e-mail: cpinto@ast. cam.ac.uk. 

2 Max-Planck-Institut fur extraterrestrische Physik, Giessenbachstrasse 1, D-85748 Garching, Germany 

3 Kavli Institute for Particle Astrophysics and Cosmology, Stanford University, 452 Lomita Mall, Stanford, CA 94305-4085, USA 

4 Department of Physics, Stanford University, 382 Via Pueblo Mall, Stanford, CA 94305-4060, USA 

5 SRON Netherlands Institute for Space Research, Sorbonnelaan 2, 3584 CA Utrecht, The Netherlands. 

6 Argelander-Institut fur Astronomie, Universitat Bonn, Auf dem Hiigel 71, 53121 Bonn, Germany. 

7 Department of Physics, University of Helsinki, FI-00014 Helsinki, Finland 

Received 5 November 2014 / Accepted 24 December 2014 


ABSTRACT 

Context. Feedback from active galactic nuclei, galactic mergers, and sloshing are thought to give rise to turbulence, which may prevent 
cooling in clusters. 

Aims. We aim to measure the turbulence in clusters of galaxies and compare the measurements to some of their structural and 
evolutionary properties. 

Methods. It is possible to measure the turbulence of the hot gas in clusters by estimating the velocity widths of their X-ray emission 
lines. The Reflection Grating Spectrometers aboard XMM -Newton are currently the only instruments provided with sufficient effective 
area and spectral resolution in this energy domain. We benefited from excellent 1.6 Ms new data provided by the Chemical Enrichment 
RGS cluster sample (CHEERS) project. 

Results. The new observations improve the quality of the archival data and allow us to place constraints for some clusters, which were 
not accessible in previous work. One-half of the sample shows upper limits on turbulence less than 500 km s _l . For several sources, 
our data are consistent with relatively strong turbulence with upper limits on the velocity widths that are larger than 1000 km s A 
The NGC 507 group of galaxies shows transonic velocities, which are most likely associated with the merging phenomena and bulk 
motions occurring in this object. Where both low- and high-ionization emission lines have good enough statistics, we find larger upper 
limits for the hot gas, which is partly due to the different spatial extents of the hot and cool gas phases. Our upper limits are larger 
than the Mach numbers required to balance cooling, suggesting that dissipation of turbulence may prevent cooling, although other 
heating processes could be dominant. The systematics associated with the spatial profile of the source continuum make this technique 
very challenging, though still powerful, for current instruments. In a forthcoming paper we will use the resonant-scattering technique 
to place lower-limits on the velocity broadening and provide further insights on turbulence. The ASTRO-H and Athena missions will 
revolutionize the velocity estimates and discriminate between different spatial regions and temperature phases. 
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1. Introduction 

Clusters of galaxies are the most massive, individual, bound 
objects in the Universe. In their gravitational potential well, 
the gas, called the intra-cluster medium (ICM), is heated to 
temperatures of 10 7 x K and, therefore, strongly emits at X- 
ray energies. The ICM is commonly thought to be in hydro¬ 
static equilibrium, but there are several factors that may af¬ 
fect the dynamical state of the gas. The feedback from active 
galactic nuclei (AGN) creates bubbles th at may drive turbu- 
lence up to about 5 00km s" 1 (see. e.g.. lBrtiggen et alj|2005l and 
Fabian et al. 2005). Sloshing of gas within the gravitational po¬ 
tential may produce similar velocities, while galactic mergers 
can give rise to even higher velocities of about 1000 km s 1 (see, 
e.g.. lLau et al.ll2009llAscasibar & Markevitchll2006l) . 

The AGN feedback is thought to offset radiative losses 
and to suppress cooling in isolated giant elliptical galax¬ 
ies and in larger syst e ms up to the richest galaxy clus¬ 
ters (see, e.g., iMcNamara & Nulsenl 120071 and iFabianl 120121) . 
Simulations and observations have confirmed that AGN feed¬ 


back may prevent cooling through the production of turbul ence 
(see, e.g., Ruszkowski et al.112004 IZhuravleva et alJ 12014 and 
iGasnari et alJ 2014 ). Other work suggests that turbulent mixing 
may also be an import ant mechanism through wh ich AGN heat 
cluster cores (see, e.g jBaneriee & SharmalooTit) . 

It is possible to measure velocity broadening on the order of 
few hundreds km s' 1 directly in the X-ray emission lines pro¬ 
duced by the hot ICM. The R eflection Grating Spectrometers 
(RGS, Iden Herder et al.l 1200 ll) aboard XMM -Newton are cur¬ 
rently the only X-ray instruments, which have enough col¬ 
lecting area and spectral resolution to enable this measure¬ 
ment. However, the spatial extent of clusters complicates the 
process due to the slitless nature of the RGS. ISanders et al.l 
(201 0) made the first measurement of cluster velocity broad¬ 
ening using the luminous cluster A 1835 at redshift 0.25. Due 
to the limited spatial extent of its b right core, an upper limit of 
274kms _1 was obtained. ISanders et alJ ( 2011 ) then constrained 
turbulent velocities for a large sample of 62 sources observed 
with XMM-AevDon/RGS, which included clusters, groups, and 
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elliptical galaxies. Ha lf of them show ve l ocit y broadening below 
700km s -1 . Recently. ISanders & Fabianl ( 2013 ) used continuum- 
subtracted emission line surface brightness profiles to account 
for the spatial broadening. This technique is affected by system¬ 
atic enorsofuptqj_50 km s - *. _ 

IWerner et al] (120091) and Ide Plaa et al] (l2012h measured tur¬ 
bulent velocities through the ratio of the Fe xvii emission lines 
at 15 and 17 A. When the velocity broadening is low, the gas is 
optically thick in the 15 A line due to resonant scattering, while 
the 17 A lines remain optically thin. The comparison of observed 
with simulated line ratios for different Mach numbers constrains 
the level of turbulence. This method is very efficient for cool 
core clusters rich in Fe xvii emission lines, but it is partly lim¬ 
ited by the systematic uncertainty (-20%) in the line ratio for an 
optically thin plasma. 

In this work, we measure the velocity broadening for the 
44 sources of the CHEmical Enrichment RGS cluster Sample 
(CHEERS), which is connected to a Very Large Program ac¬ 
cepted for XMM-Newton AO-12. We model the line spatial 
broadening using CCD images. This method has systematics due 
to the spatial profile of the continuum, which may overestimate 
the line spatial broadening, but it is still a useful technique to 
measure the level of velocity broadening when deep, high-spatial 
resolution maps are lacking. We also test an alternative method, 
which uses a variable spatial-broadening. The paper is organized 
as follows. In Sect.[2] we give a brief description of the CHEERS 
project. In Sect.[2 we present the data reduction. Our method is 
described in Sect. [4] We discuss the results in Sect. [5] and give 
our conclusions in Sect.[6] Further useful material is reported in 
AnnendixlAlto speed up the paper reading. 

2. The CHEERS project 

The current catalog includes 44 nearby, bright clusters, groups 
of galaxies, and elliptical galaxies with a value of a > 5<x detec¬ 
tion for the O viii ls-2p line at 19 A and with a well-represented 
variety of strong, weak, and non cool-core objects. This catalog 
also contains 19 new observations of 1.6 Ms in total, which are 
taken during AO-12, PI: J. de Plaa (see TableQ}. More detail on 
the sample choice is provided by another paper (de Plaa et al., in 
preparation). Among the several goals of this large project, we 
mention the following ones: 

- To understand the ICM metal enrichment by different SN 
types, (see, e.g., Mernier et al. accepted) 

- to study substructures, asymmetries and multiphaseness, 

- to study heating and cooling in cluster cores, 

- to measure turbulence (this paper), 

- to improve the cross-calibration between X-ray satellites. 

3. Data 

The data used in this paper are listed in Table|Tj In boldface, we 
show the new observations taken during AO-12. A few archival 
exposures have not been used, since they were too short. 

The XMM-Newton satellite is equipped with two types 
of X-ray detectors: The CCD-type European Photon Imaging 
Cameras (EPIC) and the Reflection Grating Spectrometers 
(RGS). The E uropean photon im agin g cameras ye M OS 1, 
MOS 2, and pn (IStriider et alJl200 H and [Turneret al.ll200ll) . The 
RGS camera consists of two similar detectors, which have both 
high effective area and spectral resolution between 6 and 38 A 
(■den Herder et alfeOOll) . The MOS cameras are aligned with the 
RGS detectors and have higher spatial resolution than the pn 
camera. We have used MOS 1 for imaging and RGS for spec¬ 
tral analysis. 



Fig. 1 . RGS extraction regions and MOS 1 stacked image of 
M 87. 


3.1. RGS and MOS 1 data reduction 

The data were reduced with the XMM-Newton Science Analysis 
System (SAS) vl3.5.0. We processed the RGS data with the SAS 
task rgsproc and the MOS 1 data with emproc to produce event 
files, spectra, and response matrices for RGS and MOS data. 

To correct for contamination from soft-proton flares, we used 
the SAS task evselect to extract light curves for MOS 1 in the 10- 
12 keV energy band, while we used the data from CCD num¬ 
ber 9 for RGS where hardly any emission from each source 
is expected. We binned the light curves in 100 s intervals. A 
Poissonian distribution was fitted to the count-rate histogram, 
and all time bins outside the 2 cr level were rejected. We built 
the good time intervals (GTI) files with the accepted time events 
for the MOS and RGS files through the SAS task tabgtigen and 
reprocessed the data again with rgsproc and emproc. The RGS 1 
total clean exposure times are quoted in Table[I] 


3.2. RGS spectra extraction 

We extracted the RGS source spectra in two alternative regions 
centered on the emission peak: a broader 3.4’ region, which in¬ 
cludes most of the RGS field of view and a narrower 0.8’ re¬ 
gion that provides the cluster cores but with high statistics. This 
was done by launching rgsproc twice by setting the xpsfincl 
mask to include 99% and 90% of point-source events inside 
the spatial source extraction mask, respectively. We have used 
the model background spectrum created by the standard RGS 
rgsproc pipeline, which is a template background file, based on 
the count rate in CCD 9. The RGS spectral extraction regions 
and the MOS 1 image of M87 are shown in Fig. [T] The spectra 
were converted to SPEXQ format through the SPEX task trafo. 
During the spectral conversion, we chose the option of sectors 
in the task trafo to create as many sectors as the different ex¬ 
posures of each source. This permits us to simultaneously fit the 
multiple RGS spectra of each source by choosing which parame¬ 
ters to either couple or unbind in the spectral models of different 
observations. 

3.3. MOS 1 spatial broadening profiles 

The RGS spectrometers are slitless, and, therefore, the spectra 
are broadened because of the spatial extent of the source in the 
dispersion direction. The effect of this spatial broadening is de¬ 
scribed by the following wavelength shift 

0.138 

A A= -A6>A, (1) 

m 

1 www.sron.nl/spex 
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Table 1. XMM-We wton /RGS observations used in this paper. 


Source 

ID ,a) 

Total clean time (ks) w 

kT (keV) <c) 

z lc) 

N H (10- 4 nr 2 ) W) 

2A0335+096 

0109870101/0201 0147800201 

120.5 

3.0 

0.0349 

30.7 

A 85 

0723802101/2201 

195.8 

6.1 

0.0556 

3.10 

A 133 

0144310101 0723801301/2001 

168.1 

3.8 

0.0569 

1.67 

A 189 

0109860101 

34.7 

1.3 

0.0320 

3.38 

A 262 

0109980101/0601 0504780101/0201 

172.6 

2.2 

0.0161 

7.15 

A 496 

0135120201/0801 0506260301/0401 

141.2 

4.1 

0.0328 

6.00 

A 1795 

0097820101 

37.8 

6.0 

0.0616 

1.24 

A 1991 

0145020101 

41.6 

2.7 

0.0586 

2.72 

A 2029 

011 1270201 0551780201/0301/0401/0501 

155.0 

8.7 

0.0767 

3.70 

A 2052 

0109920101 0401520301/0501/0601/0801 
0401520901/1101/1201/1301/1601/1701 

104.3 

3.0 

0.0348 

3.03 

A 2199 

0008030201/0301/0601 0723801101/1201 

129.7 

4.1 

0.0302 

0.909 

A 2597 

0108460201 0147330101 0723801601/1701 

163.9 

3.6 

0.0852 

2.75 

A 2626 

0083150201 0148310101 

56.4 

3.1 

0.0573 

4.59 

A 3112 

0105660101 0603050101/0201 

173.2 

4.7 

0.0750 

1.38 

A 3526 

0046340101 0406200101 

152.8 

3.7 

0.0103 

12.2 

A 3581 

0205990101 0504780301/0401 

123.8 

1.8 

0.0214 

5.32 

A 4038 

0204460101 0723800801 

82.7 

3.2 

0.0283 

1.62 

A 4059 

0109950101/0201 0723800901/1001 

208.2 

4.1 

0.0460 

1.26 

AS 1101 

0147800101 0123900101 

131.2 

3.0 

0.0580 

1.17 

AWM 7 

0135950301 0605540101 

158.7 

3.3 

0.0172 

11.9 

EXO 0422 

0300210401 

41.1 

3.0 

0.0390 

12.4 

Fornax 

0012830101 0400620101 

123.9 

1.2 

0.0046 

1.56 

HCG 62 

0112270701 0504780501 0504780601 

164.6 

1.1 

0.0140 

3.76 

Hydra-A 

0109980301 0504260101 

110.4 

3.8 

0.0538 

5.53 

M 49 

0200130101 

81.4 

1.0 

0.0044 

1.63 

M 86 

0108260201 

63.5 

0.7 

-0.0009 

2.97 

M 87 (Virgo) 

0114120101 0200920101 

129.0 

1.7 

0.0042 

2.11 

M 89 

0141570101 

29.1 

0.6 

0.0009 

2.96 

MKW 3 s 

0109930101 0723801501 

145.6 

3.5 

0.0450 

3.00 

MKW 4 

0093060101 0723800601/0701 

110.3 

1.7 

0.0200 

1.88 

NGC 507 

0723800301 

94.5 

1.3 

0.0165 

6.38 

NGC 1316 

0302780101 0502070201 

165.9 

0.6 

0.0059 

2.56 

NGC 1404 

0304940101 

29.2 

0.6 

0.0065 

1.57 

NGC 1550 

0152150101 0723800401/0501 

173.4 

1.4 

0.0123 

16.2 

NGC 3411 

0146510301 

27.1 

0.8 

0.0152 

4.55 

NGC 4261 

0056340101 0502120101 

134.9 

0.7 

0.0073 

1.86 

NGC 4325 

0108860101 

21.5 

1.0 

0.0259 

2.54 

NGC 4374 

0673310101 

91.5 

0.6 

0.0034 

3.38 

NGC 4636 

0111190101/0201/0501/0701 

102.5 

0.8 

0.0037 

2.07 

NGC 4649 

0021540201 0502160101 

129.8 

0.8 

0.0037 

2.23 

NGC 5044 

0037950101 0584680101 

127.1 

1.1 

0.0090 

6.24 

NGC 5813 

0302460101 0554680201/0301/0401 

146.8 

0.5 

0.0064 

6.24 

NGC 5846 

0021540101/0501 0723800101/0201 

194.9 

0.8 

0.0061 

5.12 

Perseus 

0085110101/0201 0305780101 

162.8 

6.8 

0.0183 

20.7 


(n) Exposure ID number. (f>l RGS net exposure time . (c) Redshifts and temperatures are adapted from| Chen et all d2007l )~an d lSnowden et all (20081 ). 


(rf) Hydrogen column density (see http://www.swift.ac.uk/analysis/nhtot/). New observations from our proposal are shown in boldface. 


where m is the spectral order and 6 is the offset angle of the 
source in arcmin (see the XMM -Newton Users Handbook). 

The MOS 1 DET Y direction is parallel to the RGS 1 disper¬ 
sion direction and can be used to correct for the spatial broaden¬ 
ing. With evselect , we extracted MOS 1 images for each expo¬ 
sure in the 0.5-1.8keV (7-25 A) energy band and their surface 
brightness profiles in the dispersion direction. We account for 
spatial broadening using the Ipro multiplicative model in SPEX, 
which convolves the RGS response with our model of the spatial 
extent of the source. We show some cumulative profiles of spa¬ 
tial broadening in Fig. [2] We have also produced stacked Fe-L 
band (10-14 A) images for each source. The central 10’ region 
contains most of the cluster emission (see Fig. IA.ll) . 


4. Spectral modeling 

Our analysis focuses on the 7 - 28 A (0.44 - 1.77 keV) first and 
second order RGS spectra. We model the spectra with SPEX 
version 2.03.03 . We scale elemental a bundances to the proto- 
Solar values of iLodders & Palmel (120091) . which are the default 
in SPEX. We adopt C-statistics and 1 cr errors throughout the pa¬ 
per, unless otherwise stated and the updated ionization balance 
calculations of lBrvans et al. d2009l) . 

Clusters of galaxies are not isothermal, and most of them 
have both hot and cool gas phases (see e.g. [Frank et aPl2013 ). 
Therefore, we have used a two-temperature thermal plasma 
model of collisional ionization emission (CIE). This model is 
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Fig. 2. MOS 1 average 7-25 A cumulative spatial profiles. 


able to fit all the spectra in our database. The cie model in 
SPEX calculates the spectrum of a plasma in collisional ioniza¬ 
tion equilibrium. The basis for this model is given by the mekal 
model, but several updates have been included (see the SPEX 
manual). Free parameters in the fits are the emission measure 
Y = n e «h dV, the temperature T, the abundances (N, O, Ne, Mg, 
and Fe), and the turbulent broadening v of the two cie models. 

We bind the parameter v and the abundances of two cie com¬ 
ponents with each other and assume that the gas phases have 
the same turbulence and abundances. This decreases the degree 
of degeneracy. This assumption is certainly not true, but some 
clusters just need one CIE component and the spectra of sev¬ 
eral clusters do not have good enough statistics in both highl¬ 
and low-ionization emission lines, which prohibits constraining 
the velocities and the abundances for both hot and cool phases. 
We attempt to constrain the turbulence in the different phases in 
Sect.ED 

The cie models are multiplied by a redshift model and a 
model for Galactic absorption, which is provided by the hot 
model in SPEX with T _= 0.5 eV and A™ r , as estimated through 
the tool of lWillingale et alj (120131) . This tool includes the contri¬ 
bution to absorption from both atomic and molecular hydrogen. 
The redshifts and column densities that have been adopted are 
shown in Table Q] To correct for spatial broadening, we have 
multiplied the spectral model by the Ipro component that re¬ 
ceives as input the surface brightness profile extracted in the 
MOS 1 images (see Sect. l3.3l and Fig. [2}. 

We do not explicitly model the cosmic X-ray background in 
the RGS spectra because any diffuse emission feature would be 
smeared out into a broad continuum-like component. 

For a few sources, including the Perseus and Virgo (M 87) 
clusters, we have added a further power-law emission compo¬ 
nent of slope ~ 2 to take the emission from the central AGN 
into account. This component is not convolved with the spatial 
profile because it is produced by a central point-like source. 

To avoid the systematic effects due to the stacking of mul¬ 
tiple observations with different pointing, we have simultane¬ 
ously fitted the individual spectra of each source extracted in the 
two regions defined in Sect. l3.2l and shown in Fig.Q] The plasma 
model is coupled between the observations. The only uncoupled 
parameters are the emission measures of the two collisional- 
ionized gas components. For each observation we adopt the spa¬ 
tial profile extracted in the MOS 1 image taken during that expo¬ 
sure. For those exposures, during which the MOS 1 detector had 
a closed filter, we have adopted an exposure-weighted average 


profile as given by the other available observations, but the 6A 
parameter in the Ipro component is left free. This factor allows 
us to shift the model lines by the same amount (in A) for each 
specific spectrum and strongly decreases the systematic effects. 
The 6A parameter is always free in our fits to account for any red- 
shift variation, which would otherwise affect the line modeling 
(see e.g. lSanders et al.ll201 ll) . 

The simultaneous modeling of multiple observations has 
been done through the use of the sectors option in SPEX (see 
also Sect. 13. 21) . The RGS 1 and 2 spectra of the same observation 
have exactly the same model and provide a single sector, while 
RGS spectra of other observations contribute additional sectors 
and have the cie normalizations uncoupled. 

4.1. Results using a fixed spatial broadening 

We have successfully applied this multi-temperature model to 
both the 3.4’ and 0.8’ RGS spectra. We show the spectral mod¬ 
eling for the 3.4’ region of the 44 sources in Figs. IA.5llA.6l and 
IA.7l in Appendix[A] We display the first-order stacked spectra to 
underline the high quality of these observations and to show the 
goodness of the modeling. 

For some sources like Fornax, M 49, M 86, NGC 4636, and 
NGC5813, the 15 and 17 A Fexvn emission lines are not well 
fitted. Precisely, the model underestimates the line peaks and 
overestimates the broadening. This may be due to the different 
spatial distribution of the gas responsible for the cool Fexvii 
emission lines and for the one producing most of the high- 
ionization Fe-L and Oviii lines. The cool gas is indeed to be 
found predominantly in the center of the clusters showing a pro¬ 
file more peaked than that one of the hotter gas. The estimated 
spatial profiles depend on the emission of the hotter gas due to 
its higher emission measure, and, therefore, they overestimate 
the spatial broadening of the 15-17 A lines. It is hard to extract a 
spatial profile for these lines because MOS 1 has a limited spec¬ 
tral resolution, and the images extrac ted in such a short ban d 
lack the necessary statistics (see e.g. ISanders & Fabianll2013h . 
In Sect. 15.11 we attempt to constrain the turbulence for lines of 
different ionization states. The 15 A/17 A line ratio is also af¬ 
fected by resonant scattering, which would require a different 
approach. We refer to a forthcoming paper on the analysis of the 
resonant scattering in the CHEERS sources. 

We skip the discussion of the abundances and the supernova 
yields because these will be treated by other papers of this series 
(de Plaa et al. in preparation and Mernier et al. submitted). 

In Fig.[3] panel), we show the upper limits on the ve¬ 
locity broadening obtained with the simultaneous fits of the 0.8’ 
7-28 A RGS spectra. We obtain upper limits for most clusters, 
while NGC 507 shows high kinematics. More detail on our re¬ 
sults for the 3.4’ and 0.8’ regions and their comparison are re¬ 
ported in Table lA. H and Fig. IA.2l (left panel). The 3.4’ limits are 
more affected by the source continuum, as clearly seen for M 87, 
AWM7, and A 4038, which makes them less reliable. 

4.2. Results using the best-fit spatial broadening 

It is known that the spatial profile of the source continuum may 
be broader than the spatial distribution of the lines. The MOS 1 
images are strongly affected by the profile of the source contin¬ 
uum and, therefore, may overestimate the spatial line broaden¬ 
ing and underestimate the residual velocity broadening. For in¬ 
stance, NGC 1316 and NGC5846 show lcr limits of 20kms _1 , 
which are not realistic (see Table lA. it . 
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Fig. 3. Left panel: Velocity 68% (red) and 90% (green) limits for the 0.8’ region with the spatial broadening determined with MOS 1 
images (see Sect. 14. 11 1. Right panel: Velocity limits obtained using the best-fit spatial broadening (see Sect l4.2l> . 


To obtain more conservative limits, we have simultaneously 
modeled the spatial and the velocity broadening. This was done 
by fitting the RGS 0.8’ spectra with a free s parameter in the 
Ipro component. This factor simply scales the width of the spa¬ 
tial broadening by a factor free to vary (see the SPEX man¬ 
ual). The free s parameter increases the degeneracy in the model 
but provides conservative upper limits on the residual velocity 
broadening, which is measured with the v parameter of the cie 
component. The new limits on the velocities are plotted in Fig. [2 
(right panel ) and quoted in the last two columns of Table lA. ll 
In Fig. IA.21 (r/g/;f panel), we compare the velocity upper limits 
estimated with the standard method (MOS 1 spatial profile with 
s = 1 in the Ipro component) with this new approach using a 
free ,v parameter. They generally agree, but the new upper limits 
on the hotter Abell clusters are systematically larger by an aver¬ 
age factor ~ 2. This confirms that some of the previous velocity 
limits were underestimated due to the broader spatial profiles. 
Therefore, we believe the new upper limits to be the most conser¬ 
vative. Interestingly, the conservative velocity limits of the hot 
clusters are generally higher than the cool galaxy groups with 
the exception of NGC 507, which is expected since the sound 
speed scales as a power of the temperature (see Sect. IO >. 


4.3. Further tests 


To estimate the contribution of the spatial broadening to the line 
widths, we have temporarily removed the convolution of the 
spectral model for the spatial profile and re-fitted the data. In 
these fits the v parameter of the cie component accounts for any 
contribution to the line broadening. The total (spatial + velocity) 
widths are also quoted in Table lA. ll 

We have also tested the co ntinuum-subtracted line s urface 
brightness profiles introduced bv ISanders & Fabianld2013i) . This 
new method consists of subtracting the surface brightness pro¬ 
files of two regions that are clearly line-dominated (core) and 
continuum-dominated (outskirts). It can be applied only to those 
objects with a narrow core where it is possible to distinguish 
between line-rich and line-poor regions. We have locally fit¬ 
ted the Ovm 19.0 A emission line of A 2597, A 3112, Hydra- 
A, Fornax (NGC 1399), and NGC 4636 and we have found a 
general agreement with the results of lSanders & Fabian s ( 201 31). 
However, our MOS 1 images have much lower spatial resolu¬ 
tion than the Chandra maps used by them, which increases the 
uncertainties that are present in this method. A thorough, exten¬ 
sive, analysis would require deep Chandra maps that are not yet 
available. 
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5. Discussion 

In this work we have analyzed the data of 44 clusters, groups of 
galaxies, and elliptical galaxies included in the CHEERS project, 
a Very Large Program that was accepted for XMM -Newton AO- 
12 (see Sect. [2]) together with complementary archival data. 

We have measured upper limits of velocity broadening for 
these objects with a me thod similar to the previo us one used by 
iBulbul et al.l (1201 2l) and lSanders & Fabianl (1201 3l) . This consists 
of fitting high-quality grating spectra by removing the spatial 
broadening through surface brightness profiles of the sources as 
provided by CCD imaging detectors. These profiles are unfor¬ 
tunately affected by the source continuum and tend to overesti¬ 
mate the line spatial broadening with a consequent s hrinking of 
the residual velocity broadening. ISanders & FabiarJ d2013l) ad¬ 
dressed this point in their Sect. 2.2 on A 3112 where they de¬ 
creased these systematic effects by using Chandra continuum- 
subtracted line spatial profiles. We have tested this method by 
using the MOS 1 observations that were taken simultaneously 
with the RGS spectra (see Fig. lA. 1 b . XMM -Newton CCDs have a 
spatial resolution lower than Chandra CCDs, which increase the 
systematic effects in the creation of continuum-subtracted maps. 
Deep Chandra observations, enabling an accurate subtraction of 
different energy bands, are missing for most sources. We have 
therefore tried to use the MOS 1 integral maps and to fit the con¬ 
tribution of the spatial broadening as an alternative method. 

5.1. Temperature dependence of the upper limits 

So far we adopted the same velocity broadening for all the emis¬ 
sion lines. For most sources it is possible to measure the velocity 
broadening of the O vm and Fe xx-to-xxiv emission lines, which 
are mainly produced by hot gas. Only a few sources have high- 
statistics Fe xvn lines produced by cool (T < 1 keV) gas. Six ob¬ 
jects exhibit both strong low- and high-ionization lines and allow 
to fit the velocity broadening of the two cie components, sepa¬ 
rately, in the full-band spectral fits. 17 sources allow to measure 
90% upper limits on turbulence for the O vm, Fe xvn, and Fe xx 
lines, by fitting the 18.0-23.0 A, 14.0-18.0 A, and 10.0-14.3 A 
rest-frame wavelength ranges, respectively. For each local fit 
we adopt an isothermal model and correct for spatial broad¬ 
ening by using additional surface-brightness profiles calculated 
through MOS 1 images extracted in the same rest-frame wave¬ 
length ranges using the same method shown in Sect. 13. 31 These 
profiles are still affected by the continuum but provide a bet¬ 
ter description of the spatial broadening in each line. In Fig. [4] 
we compare the Oviii velocity limits with those measured for 
the Fexvn and Fexx line systems. The high-ionization Fe lines 
clearly show higher upper limits, which is confirmed by the re¬ 
sults of the full-band fits: the hotter (Tl) CIE component allows 
for higher values of velocity broadening. The hotter gas is dis¬ 
tributed over a larger extent than that of the cold (T2) gas and has 
larger spatial broadening, which affects the T1-T2 results shown 
in this plot. The O vm, Fe xvn, and Fe xx lines were fitted by sub¬ 
tracting the spatial broadening extracted exactly in their energy 
band, which should partly correct this systematic effect, but it is 
difficult to estimate the systematic uncertainties due to the low 
spatial (and spectral) resolution of the CCD data. On some ex¬ 
tent, the hotter phase may still have larger turbulence. For clarity, 
we also tabulate these line-band fits in Table lA.2l We also note 
that the the velocity limits of low- and high-ionization iron lines 
fall at opposite sides of the Fe-0 vm 1:1 line, which means that 
the metallicity distribution in the sources should not affect our 
broad-band, multi-ion, limits shown in Fig. [3] 



Fig. 4. 90% upper limits on velocity broadening obtained in the 
0.8’ region for the Ovm lines compared with those measured 
for high-ionization Fe xx (open red triangles) and for the low- 
ionization Fexvn (filled black triangles) line systems. For six 
sources we could also measure the limits for the hot (open green 
circles) and the cool (open magenta boxes) CIE components (see 
Scci.FTl). 


5.2. Turbulence 

In Fig. [3] we show the velocity broadening of the RGS spectra 
extracted in the 0.8’ core region. We find upper limits to the ve¬ 
locity broadening with the possible exception of NGC 507. They 
generally range between 200 and 600km s 1 . For several objects 
like A85, A 133, M49, and most NGC elliptical we found ve¬ 
locity levels below 500 km s _1 , which would suggest low turbu¬ 
lence. The broader 3.4’ region is more affected by spatial broad¬ 
ening as shown by the higher upper limits, but non-detection, 
of A 4038, AWM7, and M 87 (see TablelAJland Fit'.fACl For 
these sources it is difficult to constrain the velocity broadening 
because their large extent smears out the emission lines. 

To understand how much energy can be stored in turbulence, 
we compare our upper limits with the sound speeds and the 
temperatures of dominant cie component in these objects. The 
sound speed is given by cs = y/yk T/pm v , where y is the adi¬ 
abatic index, which is 5/3 for ideal monoatomic gas, T is the 
RGS temperature, p = 0.6 is the mean particle mass, and m p 
is proton mass. The ratio between turbulent and thermal energy 
is fiturb/stheim = y/2 M 2 , where M = v tur b/<A is the Mach num¬ 
ber (see also IWerner et al.l 1 2009 ). In Fig. lA.21 we compare our 
2cr upper limits on the velocities in the central 0.8’ region with 
the sound speed and some fractions of turbulent energy. We also 
show the more conservative velocity upper limits that were mea¬ 
sured with a variable spatial broadening. At least for half the 
sample, our 90% upper limits are below the sound speed in the 
system. In about ten objects the turbulence contains less than the 
4 0% of the thermal energ y. This is similar to the previous results 
of ISanders & Fabianl (120131) . Apparently, the hotter objects allow 
for higher velocities. 

We note that the spectral extraction region had a fixed angle, 
and, therefore, the actual physical scale - where we estimated 
the velocity broadening - depends on the source distance. In 
Fig. lA.41 (left panel), we show the Mach numbers for the 90% 
conservative upper limits as a function of the temperature, and 
we compare the average upper limits on Mach number calculated 


























Ciro Pinto et al.: Chemical Enrichment RGS cluster sample (CHEERS): Constraints on turbulence 


7 


within different ranges of physical scales. There is no significant 
trend with the temperature, but the average upper limit on the 
Mach number is lower for narrower physical scales. Assuming 
a Kolmogorov spectrum for the turbulence in these objects, the 
root-mean-square velocity scale depends on the l/3rd power of 
the physical length. Therefore, we scaled the upper limits by 
(■sc/sCmjn) 1 ^, where sc m i n is the minimum physical scale per arc- 
sec ~ 0.07 kpc/1” of NGC 4636, the nearest object in our sample. 
In other words, we divided our upper limits by the relative physi¬ 
cal scale per arcsec relative to NGC 4636, which is equivalent to 
normalizing by the ratio between the size of the spectral extrac¬ 
tion region of each cluster and that one of NGC 4636. The scaled 
upper limits on the Mach numbers are tabulated in Table lA.2l and 
plotted in F\s \AA\ (rieht panel). They are randomly distributed 
around Ma ~ 0.8 and do not depend any more on the physical 
scale. We coded the point-size and the colors with the values of 
r 5 oo and K {) taken from the literature. The 7*500 is the radius within 
which the mean over-density of the cluster is 500 times the criti¬ 
cal density at the cluster redshift, and Kq is the value of the cen¬ 
tral entropy in the same cluster. All the adopted values and their 
references are reported in Table lA.2l We do not find any signifi¬ 
cant relation between the upper limits on the Mach number and 
these physical properties, possibly due to the limited sample. 

To understand whether dissipation of turbulence may pre¬ 
vent cooling in our sample, we computed the Mach number that 
is required to balance the heating and cooling, according to the 
following equation: 


Ma 


REQ 


) !5 ( v /3 

\10- 2 cm- 3 / 


10 3 kms 


U 

lOkpc / 


1/3 


( 2 ) 


where n e is the density at the cavity location, c s the sound speed 
that we have estimated through the RGS temperature, and / the 
chara cteristic eddy size, w hich we take as the average cavity size 
(see lZhuravleva et~ai]l2014l) . The Mach numbers required to bal¬ 
ance cooli ng are tabulated in Table lA.2l Most cavity sizes were 
taken from lPanagoulia et al.l d2014bh . For clusters with multiple 
cavities, we used an average size. For the 19 sources outside of 
their sample, we used their r — T relation to determine the cavity 
size. Most densities were taken from the ACCEPT catalog. 

In Fig. [3 we compare the ratios between the conservative 
upper limits of the scaled Mach numbers assuming Kolmogorov 
turbulence, and those that are required to balance cooling with 
the RGS temperatures. For most sources, our upper limits are 
larger than the balanced Mach numbers, which means that dis¬ 
sipation of turbulence can provide enough heat to prevent the 
cooling of the gas in the cores. 

It is difficult to know which is the main mechanism that pro¬ 
duces turbulence in these objects. Our scaled upper limits are 
mostly below 500 km s -1 , which can be produced by bubbles 
inflated by past AGN activity (see, e.g.. iBriiggen et akll2005l) . 
For some objects, our upper limits are consistent with veloci¬ 
ties up to 1000 k m s 1 , which would correspond to Mach num¬ 
bers larger than one. For NGC 507, we detect transonicmotions 
presumably due to merging (see, e.g., iAscasibar & Markevitchi 
12 006 ). In a forthcoming paper, we will analyze the resonant scat¬ 
tering of the Fe xvn lines exhibited by half of our sample to place 
lower limits on turbulent broadening and provide more insights 
on its origin and its role in preventing cooling. 


5.3. Comparison with previous results 

Our velo city limits broadly agree w ith the previous results ob¬ 
tained by ISanders & Fabian! ( 2013 ) using a similar method and 



Fig. 5. Ratios between the 90% conservative upper limits on the 
Mach number (velocity / sound speed) that are scaled by the 
l/3rd power of the spatial scalel assuming Kolmogorov turbu¬ 
lence (see Sect. E3, and the Mach number, which is required 
to make a heating-cooling balance (see Eq.[2]). The point size 
provides the 7 * 500 , and the color is coded according to the central 
entropy, Kq, in units of keV cm 2 . 


by oth er authors, who use the measurements of re sonant scat¬ 
tering (IWerner et al.ll2009l and Ide Plaa et afll2012l) . In particu¬ 
lar, our limits for M49 (also known as NGC 4472), NGC 4636, 
and NGC 5813 agree with the 1001<m s 1 upper limit ob¬ 
tained by IWerner et al.l ( 20091 ). We also found upper limits of 
a few lOOskms -1 for A 3112, which is similar to the results 
of iBulbul et al.l (120121) . However, we measured higher limits 
with a variable spatial broad ening that agree with c ontinuum- 
subtracted profiles method of Sanders & F abianl (120131) . 

Recently, IZhuravleva et al.l (12014l) used the surface bright¬ 
ness fluctuations in the Chandra images of the Perseus and Virgo 
clusters to derive turbulent velocities in the range 70-210km s 1 
for Perseus and 43-140 km s _1 for Virgo, where the smaller val¬ 
ues refer to the central 1.5’ region. Our upper limits in the cores 
of the clusters are consistent with their values, especially when 
normalized by the physical scale factor 1. 570 . 4 '. They show that 
these turbulent motions should dissipate enough energy to offset 
the cooling of the central ICM in these clusters. For ten objects, 
the scaled Mach number can be transonic, and a major fraction 
of energy can be stored in turbulence, which could significantly 
heat the gas th rough dissipat i on (se e, e.g., iRuszkowski et al.l 
2004). Recently. iGasnari et al] (120141) noted that even if the tur¬ 
bulence in the hot gas is subsonic, i t may be transonic in the 
cooler gas phases. IZhuravleva et al. ( 2014 ) reported that dissi¬ 
pation of turbulence may balance cooling even under subsonic 
regime. Our upper limits on Mach number are larger than the 
values necessary to balance cooling and are consistent with this 
scenario. However, it is possible that other processes are dom¬ 
inant, such as turbulent mixing (see, e.g., iBaneriee & Sharmal 
120141) . 

The NGC 507 group exhibits velocities larger than 
10001<m s~ 1 in both the 0.8’ and 3.4’ regions, corresponding to 
a scaled Mach number Ma = 4.2 ± 1.7 (lcr). The 15 A Fexvn 
line is stronger than the one at 17 A, which would suggest low 
resonant scattering (see Fig. lA.7l > and, therefore, high kinemat¬ 
ics in the galaxy group. This object is known to have a dis¬ 
turbed shape and to host radio lobes presumably in a transonic 
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Fig. 6. Line-of-sight velocity versus projected distance from the 
central cD galaxy for the member galaxies of NGC 507 group. 
Optical spectroscopic redshifts are taken from Zhang et al. 

m . 


expansion/inflation (iKraft et al.l [2004}) . However, our high val¬ 
ues suggest the presence of bulk motions. In Fig. [6] we show the 
velocities of the galaxies in the NGC 507 group as taken from 
IZhang et~akl (1201 ll ). They are not necessary linked to that of the 
ICM, but there are high kinematics and hints of infalling clumps, 
which indicate a substructure extended toward the observer. In 
this group, the galaxy velocities generally double those observed 
in NGC 4636, where we measure lower velocity broadening (see 
also the different line widths in Fig. lA.7l >. 

5.4. Toward ASTRO-H 

The RGS gratings aboard XMM-Newton are currently the only 
instruments that can measure 100s km s 1 velocities in X-ray 
spectra of extended sources like clusters of galaxies. However, 
they are slitless spectrometers and, therefore, affected by spa¬ 
tial broadening. We have partly solved this issue by using line 
surface brightness profiles, but there are still systematic un¬ 
certainties larger than 100 kms -1 . Our models provide an im¬ 
portant workbench once the new ASTRO-H X-ray satellite 
iTakahashi et al.l l2 0ld ) is launched. The spectra, as provided by 
its microcalorimeter (SXS), do not suffer from spatial broaden¬ 
ing as for the RGS and will revolutionize the method. Moreover, 
its constant spectral resolution in terms of energy increases 
the sensitivity at high energies, which allows us to use higher- 
ionization lines up to 6-7 keV (Fe-K line complex) necessary to 
constrain the turbulence in hotter gas phases. The position of the 
lines unveil evidence of bulk motions. 

In Fig. [7] (left panel), we compare the effective area of the 
ASTRO-H SXS with that of the first order RGS 1 and 2. 
ASTRO-H provides clearly better results than the sum of RGS 1 
and 2 below 14 A (above 1 keV). The RGS has still a better 
spectral resolution than the SXS in the wavelength range that in¬ 
cludes the Fe xvn lines of the cool gas, but the absence of spuri¬ 
ous line-broadening in the SXS makes it a great alternative tool. 
We have simulated a lOOks exposure with the ASTRO-H SXS 
for four interesting objects in our catalog: Perseus (500 km s '), 


NGC5846 (lOkms- 1 ), NGC4636 (100 kins’ 1 ), and NGC507 
(1000 km s , see Fig. [7 ] right panel). We have used the model 
fitted for the full (-1.7’,+1.7’) RGS spectra as a template, which 
are shown in Fig. IA.7l because this extraction region is compa¬ 
rable to the 3.05’ x3.05’ field-of-view of the microcalorimeter. 
The spatial broadening was excluded from the model. The sim¬ 
ulated SXS spectra are characterized by a richness of resolved 
emission lines, which provides velocity measurements with an 
accuracy of 50 km s 1 or better. The line widths clearly increase 
throughoutNGC 5846, NGC 4636, and NGC 507. The hotter gas 
present in the Perseus cluster produces strong higher-ionization 
lines above 1 keV, which constrain the turbulence in different 
(Fe-L and Fe-K) gas phases. 

We also note that the 1 ’ spatial resolution of ASTRO-H pro¬ 
vides, for the first time the means for a spatially-resolved high- 
resolution spectral analysis and the measurements of turbulence 
in different regions of the clusters. The ATHENA X-ray obser¬ 
vatory that is to be launched by the late 2020s will further rev¬ 
olutionize our measurements due to its combined high spectral 
(2.5 eV) and spatial (< 5”) resolution. 


6. Conclusion 

We have presented a set of upper limits and measurements of the 
velocity widths for the soft X-ray emitting gas of a sample of 
clusters, groups of galaxies, and elliptical galaxies included in 
the CHEERS project. We have subtracted the instrumental spa¬ 
tial broadening through the use of surface brightness profiles ex¬ 
tracted in the MOS 1 images. 

For most sources, we obtain upper limits ranging within 200- 
600 km s’ 1 , where the turbulence may originate in AGN feed¬ 
back or sloshing of the ICM. However, for some sources, such 
as NGC 507, we find upper limits of 1000 km s~ 1 or larger, sug¬ 
gesting other origins, such as mergers and bulk motions. The 
measurements depend on the angular scale and the temperature. 
For a small sample producing strong high- and low-ionization 
lines, we measured significantly broader upper limits for the hot 
gas phase, which may be partly due to its larger spatial extent as 
compared to the cool phase. When we normalize the Mach num¬ 
bers for the physical scale, assuming Kolmogorov turbulence, 
we constrain upper limits ranging within 0.3 < Ma <1.5. These 
values are above the Mach numbers necessary to balance cool¬ 
ing, which means that the dissipation of turbulence may be the 
dominant mechanism to heat the gas and quench cooling flows. 
However, it is possible that additional processes are heating the 
ICM. In a forthcoming paper, we will use the resonant-scattering 
technique to place some lower limits on the velocities in half of 
our sample, which are characterized by strong Fe xvn emission 
lines, and make an extensive use of higher-resolution Chandra 
maps. This will provide alternative measurements and further in¬ 
sights on the origin and the role of turbulence in clusters, groups 
of galaxies, and elliptical galaxies. 

The current techniques are partly limited by systematics as¬ 
sociated with the spatial line broadening, but the future ASTRO- 
H and ATHENA missions will dramatically improve the method 
due to their high broadband spectral resolution and to the ab¬ 
sence of issues concerning the spurious spatial broadening. They 
will also enable spatially-resolved turbulence measurements. 
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Fig. 7. Left panel : RGS 1 and 2 first order and ASTRO-H SXS effective area. Right panel'. ASTRO-H SXS lOOks simulations for 
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Appendix A: Maps, spectra, and velocity limits 

We have placed the MOS 1 images, the RGS spectra, and the 
velocity limits in this section to unburden the paper reading. 


The MOS 1 images (Fig. IA.Il l are obtained by stacking the 
Fe-L (10-14 A) band images extracted in each exposure (see 
Sect. 13. 3k We zoomed on the central 10’ region. 

In Table lA.ll we quote all our velocity results. In Fig. IA.2l 
(left panel), we compare the 2 cr upper limits on the velocity 
broadening as measured for the 3.4’ and 0.8’ regions by sub¬ 
tracting the MOS 1 spatial profiles. In Fig. IA.2l (rightpanel), we 
compare the 0.8’ 2 cr upper limits estimated by subtracting the 
adopted and the best-fit, spatial-line-broadening (see Sect. 113- 
In Table lA. ll we also show the total line widths as result of spa¬ 
tial plus Doppler broadening with their 68% uncertainties. These 
total widths are clearly dominated by the spatial broadening. 

In Tablef/Ol we report the values of r $oo and Kq adopted 
with with the physical scales, the RGS temperatures estimated 
with an isothermal model, the conservative upper limits on the 
velocities and the Mach numbers, the limits separately mea¬ 
sured for the O vm, Fe xvn, and Fe xx-to-xxiv emission lines, 
and finally the velocity limits for the two CIE components 
where it was possible to fit them separately. The separate Fe xvn 
and 2-T fits were accessible only for a very limited sample of 
sources with both strong high- and low-ionization Fe lines (see 
Sect. 15. Ik In Fig. IA.3l we compare the velocity upper limits es¬ 
timated with both the standard and the conservative methods 
with the average RGS temperature measured with an isother¬ 
mal model, the sound speed, and the fractions of thermal en¬ 
ergy stored in turbulence. In Fig. lA.4l (/e/f panel), we show the 
Mach number as a function of the temperature with the point- 
size and, it is color coded according to the physical scale and the 
central entropy, respectively. The lines show the average Mach 
number calculated within particular ranges of physical scales. In 
Fig. IA.4l (/e/7 panel), we show the Mach number scaled by the 
l/3rd power of the physical scale, assuming Kolmogorov turbu¬ 
lence (see also Sect. E3. 

In Figs. IA.5l IA.6I and IA.7I we show the RGS spectra ex¬ 
tracted in the 3.4’ cross-dispersion region (see also Sect. 13. 2k 
The spectra were combined with the rgscombine task for plot¬ 
ting purposes. We have adapted the spectral model from the 
bestfit obtained for the parallel modeling of the individual expo¬ 
sures. Most spectra were fitted with a multi-temperature, two-cie 
model (see Sect. S. A few of them require just a single cie com¬ 
ponent. The spectral modeling involved the 7 - 28 A band, but 
we focus on the shorter 10 - 21 A band containing the Fe, Ne, 
and O lines. 
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NGC 5846 + 


5846 Perseus 


Fig.A.l. MOS 
Sect. 153V 


1 stacked Fe-L band images: Central 10’ x 10’ region. The starred clusters are part of our new campaign (see also 




Fig. A.2. Left panel : Velocity broadening at 2 cr upper limits for the (-1.7’,+ 1.7’) and the (-0.4’,+0.4’) regions at comparison. The 
spatial broadening was removed through the MOS 1 surface brightness profiles. Right panel. (-0.4’,+0.4’) velocity 2 cr upper limits 
compared with those estimated in the same region but with the variable best-fit, spatial broadening (scale parameter, s, is free in the 
Ipro component, see Sect. l4.2l and Table lA. it . 
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Fig.A.3. 90% upper limits on velocity broadening obtained in the 0.8’ region versus RGS temperature (the red arrows provide the 
conservative limits measured with the best-fit spatial broadening, see Sect. 14.21) . The sound speed and the fractions of thermal energy 
in turbulence are shown. 
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Fig. A.4. Left panel. 90% conservative upper limits on Mach number versus temperature. Point size refers to the physical scale and 
the color is coded according to the central entropy, Ko, in units of keV cm" 2 . The lines show the upper limits on Mach number 
averaged in the four ranges of physical scales to underline its dependence on the source distance. Right panel : Same as before, but 
here the limits on Mach number are scaled by the 1/3rd power of the physical scale for a comparison between sources at different 
redshift (assuming Kolmogorov turbulence, see Sect. 1531 . 
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Table A.l. Velocity broadening upper limits and total line widths. 



3.4’ 

RGS extraction region 

0.8’ RGS extraction region 

0.8’ free spatial broadening 

Source 

l^lo) 

2a ^ 

FWHM (c) 

1 cr ( “) 

2a ^ 

FWHM (C > 

1 o-»> 

2a- w 

2A0335+096 

80 

240 

2920 ± 140 

150 

400 

2240 ± 140 

570 

880 

A 85 

290 

440 

2100 ± 500 

220 

340 

720 ± 160 

580 

910 

A 133 

280 

440 

1610 ± 220 

230 

380 

1170 ± 140 

980 

1050 

A 189 

530 

840 

1980 ± 320 

700 

910 

1350 ± 190 

600 

950 

A 262 

180 

370 

4570 ±410 

180 

310 

2340 ± 450 

170 

300 

A 496 

310 

540 

3150 ±220 

270 

450 

1970 ± 170 

1250 

1520 

A 1795 

610 

990 

2590 ± 320 

540 

810 

1687 ± 320 

2100 

2300 

A 1991 

270 

500 

830 ± 250 

250 

460 

480 ± 230 

300 

520 

A 2029 

330 

560 

1470 ± 700 

410 

690 

720 ± 260 

330 

500 

A 2052 

190 

420 

2800 ± 500 

180 

430 

2289 ± 220 

570 

910 

A 2199 

380 

620 

3800 ± 300 

240 

420 

2540 ± 230 

330 

560 

A 2597 

480 

720 

1730 ± 270 

300 

510 

1130 ± 170 

1200 

1300 

A 2626 

900 

1500 

3300 ± 900 

850 

1250 

2440 ± 610 

1600 

2300 

A 3112 

420 

650 

2000 ± 170 

270 

450 

1120 ± 130 

1240 

1390 

A 3526 

260 

360 

3930 ± 110 

190 

260 

2470 ± 160 

290 

400 

A 3581 

610 

980 

4390 ± 220 

410 

670 

2980 ± 260 

700 

960 

A 4038 

5180 

6480 

8100 ± 800 

930 

1540 

4700 ± 750 

2600 

3200 

A 4059 

420 

680 

2620 ± 250 

360 

600 

1860 ± 180 

1490 

1730 

AS 1101 

340 

550 

2650 ± 200 

270 

450 

1790 ± 180 

600 

910 

AWM 7 

8000 

10000 

12600 ± 700 

1000 

1430 

10500 ± 700 

150 

600 

EXO 0422 

480 

770 

1520 ± 280 

400 

680 

1020 ± 230 

1300 

1400 

Fornax 

160 

260 

2500 ± 500 

160 

230 

820 ± 100 

170 

260 

HCG 62 

150 

290 

2180 ± 190 

140 

270 

1420 ± 90 

350 

560 

Hydra-A 

360 

730 

3990 ± 400 

550 

920 

2420 ± 370 

460 

1050 

M 49 

190 

310 

3180 ± 190 

150 

240 

1350 ± 100 

220 

320 

M 86 

650 

1100 

7500 ± 400 

400 

690 

5960 ± 490 

390 

660 

M 87 

3800 

4100 

5290 ± 150 

100 

210 

3850 ± 150 

170 

470 

M 89 

740 

910 

1700 ± 180 

710 

900 

1300 ± 190 

460 

670 

MKW 3 s 

350 

600 

4010 ± 570 

360 

590 

2170 ± 250 

430 

680 

MKW 4 

310 

550 

3620 ± 330 

360 

580 

2040 ± 140 

380 

620 

NGC 507 

2500 ± 1000 

2500 ± 1500 

6100 ± 900 

1500 + 800 

iJW — 1300 

2800 

4600 ± 500 

3000 ± 1200 

5000 

NGC 1316 

20 

210 

1770 ± 120 

no 

330 

1160 ±90 

90 

240 

NGC 1404 

250 

390 

1930 ± 130 

520 

680 

1470 ± 130 

540 

650 

NGC 1550 

860 

1320 

7500 ± 700 

890 

1100 

3750 ± 250 

1130 

1250 

NGC 3411 

290 

710 

3830 ± 360 

610 

1030 

2230 ± 310 

370 

720 

NGC 4261 

160 

300 

760 ± 120 

140 

220 

470 ± 350 

160 

230 

NGC 4325 

300 

520 

1570 ± 240 

310 

490 

1060 ± 180 

440 

680 

NGC 4374 

120 

250 

1770 ± 140 

180 

310 

1400 ± 490 

160 

280 

NGC 4636 

140 

190 

2410 ± 70 

130 

180 

1660 ± 60 

190 

300 

NGC 4649 

210 

320 

1520 ± 80 

130 

250 

1050 ± 70 

280 

390 

NGC 5044 

80 

220 

4170 ± 180 

100 

240 

3100 ±200 

140 

400 

NGC 5813 

100 

140 

4080 ± 80 

120 

160 

2980 ± 120 

110 

160 

NGC 5846 

40 

80 

2100 ± 60 

20 

130 

1670 ±50 

100 

190 

Perseus 

320 

460 

4950 ± 140 

620 

880 

4381 ±390 

1000 

1600 


Velocity lcr and 2cr upper limits for the RGS spectra extracted within the 3.4’ and 0.8’ regions centered on the source emission peak. 

For NGC 507, we quote the best-fit velocities. The spatial line broadening was subtracted through the MOS 1 maps (see Figs.lTland lA^I /e/i panel). 
(f,) Velocity limits for the 0.8’ region obtained by subtracting the best-fit spatial line broadening (see Sect. l4.2l and Fig. lA,2l n'efa panel). 

(c) Total (spatial + velocity) line widths with lcr uncertainties (see Sect. l4.3t . 
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Table A.2. Additional results and physical properties. 


Source 

d(Mpc) 

Physical properties 
kpc/1” r 50 o K 0 T 

Cs 

Ma RE Q 

RGS-band 2cr limits (h> 
v Vsc Ma Ma sc 

Emission-line 2 a limits w 

O vm Fe xvn Fe xx+ 

2-T 2 cr limits (i) 

Vl v 2 

2A0335+096 

149 

0.72 

1.05 

7.1 

1.6 

640 

0.27 

880 

400 

1.37 

0.63 

580 

- 

530 

- 

- 

A 85 

238 

1.15 

1.21 

12.5 

2.3 

780 

0.28 

910 

360 

1.16 

0.46 

480 

- 

570 

- 

- 

A 133 

243 

1.18 

0.94 

17.3 

4.5 

1090 

0.13 

1050 

410 

0.96 

0.37 

460 

- 

603 

- 

- 

A 189 

137 

0.66 

0.50 

4.0 

0.9 

500 

0.18 

950 

450 

1.91 

0.90 

650 

- 

3010 

- 

- 

A 262 

69 

0.34 

0.74 

10.6 

1.7 

670 

0.17 

300 

180 

0.45 

0.27 

510 

- 

790 

- 

- 

A 496 

140 

0.68 

1.00 

8.9 

2.8 

870 

0.18 

1520 

710 

1.75 

0.82 

480 

- 

710 

- 

- 

A 1795 

264 

1.28 

1.22 

19.0 

1.8 

690 

0.36 

2300 

870 

3.35 

1.27 

930 

- 

1750 

- 

- 

A 1991 

251 

1.22 

0.82 

1.5 

2.7 

840 

0.21 

520 

200 

0.62 

0.24 

630 

- 

900 

- 

- 

A 2029 

328 

1.59 

1.33 

10.5 

2.6 

840 

0.23 

500 

180 

0.60 

0.21 

810 

- 

2310 

- 

- 

A 2052 

149 

0.72 

0.95 

9.5 

3.2 

920 

0.21 

910 

420 

0.99 

0.45 

600 

- 

690 

- 

- 

A 2199 

129 

0.63 

1.00 

13.3 

1.3 

590 

0.28 

560 

270 

0.95 

0.46 

540 

- 

1010 

- 

- 

A 2597 

365 

1.77 

1.11 

10.6 

2.6 

840 

0.23 

1300 

440 

1.55 

0.53 

500 

- 

990 

- 

- 

A 2626 

245 

1.19 

0.84 

23.2 

1.3 

600 

0.23 

2300 

890 

3.85 

1.49 

- 

- 

- 

- 

- 

A 3112 

315 

1.53 

1.13 

11.4 

1.4 

600 

0.25 

1390 

500 

2.30 

0.82 

500 

- 

570 

- 

- 

A 3526 

44 

0.21 

0.83 

2.3 

2.3 

790 

0.20 

400 

280 

0.51 

0.35 

390 

230 

730 

530 

240 

A 3581 

91 

0.44 

0.72 

9.5 

2.7 

850 

0.14 

960 

520 

1.13 

0.61 

1150 

- 

1060 

- 

- 

A 4038 

128 

0.62 

0.89 

37.9 

2.3 

790 

0.18 

3200 

1540 

4.06 

1.96 

1540 

- 

5290 

- 

- 

A 4059 

215 

1.04 

0.96 

7.1 

3.2 

920 

0.19 

1730 

700 

1.88 

0.76 

1090 

- 

1120 

- 

- 

AS 1101 

248 

1.20 

0.98 

10.4 

2.4 

800 

0.23 

910 

350 

1.15 

0.44 

530 

- 

690 

- 

- 

AWM 7 

73 

0.35 

0.86 

8.4 

1.9 

700 

0.22 

600 

350 

0.85 

0.50 

1650 

- 

- 

- 

- 

EXO 0422 

167 

0.81 

0.89 

13.8 

2.4 

800 

0.18 

1400 

620 

1.75 

0.77 

1140 

- 

1170 

- 

- 

Fornax 

19 

0.09 

0.40 

2.6 

1.0 

520 

0.16 

260 

240 

0.50 

0.46 

390 

300 

520 

- 

- 

HCG 62 

60 

0.29 

0.46 

3.4 

0.8 

480 

0.23 

560 

350 

1.18 

0.73 

910 

320 

1350 

2370 

260 

Hydra-A 

222 

1.08 

1.07 

13.3 

3.4 

950 

0.21 

1050 

420 

1.10 

0.44 

700 

- 

1050 

- 

- 

M 49 

17.1 

0.08 

0.53 

0.9 

0.9 

490 

0.17 

320 

300 

0.65 

0.62 

520 

340 

940 

1050 

240 

M 86 

16 

0.08 

0.49 

8.0 

0.8 

450 

0.22 

660 

640 

1.46 

1.41 

- 

- 

- 

- 

- 

M 87 

16.4 

0.08 

0.75 

3.5 

1.4 

600 

0.25 

470 

450 

0.78 

0.75 

400 

- 

430 

- 

- 

M 89 

15.3 

0.07 

0.44 

3.0 

0.6 

400 

0.20 

670 

660 

1.67 

1.63 

1330 

960 

1430 

- 

- 

MKW 3 s 

192 

0.93 

0.95 

23.9 

2.7 

840 

0.16 

680 

290 

0.81 

0.34 

720 

- 

1310 

- 

- 

MKW 4 

80 

0.39 

0.62 

6.9 

1.5 

630 

0.31 

620 

350 

0.99 

0.56 

660 

- 

700 

- 

- 

NGC 507 

65.5 

0.32 

0.60 

1.0 

0.7 

420 

0.22 

5000 

3010 

11.87 

7.15 

4210 

1310 

4560 

6540 

1680 

NGC 1316 

19 

0.09 

0.46 

1.0 

0.7 

420 

0.20 

240 

220 

0.57 

0.52 

410 

250 

4780 

- 

- 

NGC 1404 

20 

0.10 

0.61 

5.5 

1.2 

560 

0.23 

650 

580 

1.16 

1.04 

1020 

690 

1180 

- 

- 

NGC 1550 

53.6 

0.26 

0.62 

6.6 

0.9 

500 

0.22 

1250 

800 

2.49 

1.60 

860 

- 

1010 

- 

- 

NGC 3411 

66.3 

0.32 

0.47 

3.7 

0.7 

440 

0.21 

720 

430 

1.65 

0.99 

1790 

1520 

1570 

- 

- 

NGC 4261 

29.4 

0.14 

0.45 

3.0 

0.9 

500 

0.19 

230 

180 

0.47 

0.37 

650 

290 

620 

- 

- 

NGC 4325 

108 

0.52 

0.58 

4.7 

0.7 

420 

0.29 

680 

350 

1.63 

0.83 

1690 

540 

1870 

- 

- 

NGC 4374 

18.4 

0.09 

0.46 

8.0 

0.7 

420 

0.19 

280 

260 

0.66 

0.61 

510 

550 

900 

- 

- 

NGC 4636 

14.3 

0.07 

0.35 

1.4 

0.9 

480 

0.13 

300 

300 

0.62 

0.62 

310 

290 

1410 

- 

- 

NGC 4649 

17.3 

0.08 

0.53 

4.6 

0.9 

490 

0.19 

390 

370 

0.80 

0.75 

370 

580 

590 

- 

- 

NGC 5044 

38.9 

0.19 

0.56 

2.3 

1.1 

540 

0.24 

400 

290 

0.74 

0.53 

500 

350 

1770 

790 

280 

NGC 5813 

29.7 

0.14 

0.44 

1.4 

0.7 

430 

0.24 

160 

130 

0.37 

0.29 

470 

170 

1690 

- 

- 

NGC 5846 

26.3 

0.13 

0.36 

1.8 

0.7 

430 

0.39 

190 

160 

0.44 

0.36 

320 

190 

1470 

300 

180 

Perseus 

73.6 

0.36 

1.29 

19.4 

2.8 

860 

0.28 

1600 

926 

1.86 

1.08 

1160 

- 

1670 

- 

- 


?3nn (in Mpc units) an d K 0 (in keV cm 2 units) were ta ken from the ACCEPT catalog dCavagnolo et al.ll2 009l) with the exceptions of M 86 
dFinoguenov et all 12004) . NGC1316 dTashiro et all 12006). NGC4649, and NGC4261 dWemer et all 120 2l). A 1 89. Fornax, HCG 62, M49, 
NGC1550, NGC3411, NGC4325, NGC4636, NGC5044, NGC5813, and NGC5846 iPanagoulia et al.ll20l4al). NGC4374 and NGC 1404 


profiles were calculated following the method of IWerner et all J2012i) . The temperatures (in keV units) are estimated with the 0.8’ region RGS 
best-fit isothermal model, and they are used to compute the sound speed, c s , which is reported in units of km s^ 1 (see Sect. 1521 
Ma RE Q is the Mach number which is required to make a heating-cooling balance. It is calculated with Eq. 101 in Sect. 15.21 

90% limits on velocity and Mach number for the 0.8’ region, obtained by subtracting the best-fit spatial line broadening, and 90% limits 
normalized by the physical scale, assuming Kolmogorov turbulence (see Sects. l4.2l and l5.21 Figs.[3land lA.31 and last column in Table lA.ll l. 

(c) 90% velocity limits separately measured for the relevant emission lines with an isothermal cluster model (see Sect. ED . 

(d) 90% velocity limits separately measured for the two CIE components with the multi-temperature cluster model (see Sect. 1511 






















Fig. A.5. RGS spectral fits for the (-1.7’,+ 1.7’) region with the 7 - 25 A spatial broadening profile (Part I). 
For displaying purposes, the spectra were combined using the XMM-SAS task rgscombine. 
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Fig. A.6. RGS spectral fits for the (-1.7’,+ 1.7’) region with the 7 - 25 A spatial broadening profile (Part II). 
For displaying purposes, the spectra were combined using the XMM-SAS task rgscombine. 
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